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ABSTRACT 

The formation of gas giant planets is assumed to be preceded by the emer- 
gence of sohd cores in the conventional sequential-accretion paradigm. This hy- 
pothesis implies that the presence of earth-like planets can be inferred from the 
detection of gas giants. A similar prediction cannot be made with the gravita- 
tional instability (hereafter GI) model which assumes that gas giants (hereafter 
giants) formed from the collapse of gas fragments analogous to their host stars. 
We propose an observational test for the determination of the dominant planet- 
formation channel. Based on the sequential-accretion (hereafter SA) model, we 
identify several potential avenues which may lead to the prolific formation of a 
population of close-in earth- mass (M^) planets (hereafter close-in earths) around 
stars with 1) short-period or 2) solitary eccentric giants and 3) systems which 
contain intermediate-period resonant giants. In contrast, these close-in earths 
are not expected to form in systems where giants originated rapidly through GI. 
As a specific example, we suggest that the SA processes led to the formation of 
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the 7.5 M® planet around GJ 876 and predict that it may have an atmosphere 
and envelope rich in O2 and liquid water. Assessments of the ubiquity of these 
planets will lead to 1) the detection of the first habitable terrestrial planets, 2) 
the verification of the dominant mode of planet formation, 3) an estimate of the 
fraction of earth-harboring stars, and 4) modification of bio-marker signatures. 

Subject headings: planetary systems: formation - planetary systems: protoplan- 
etary disks - planets and satellites: formation - solar system: formation - stars: 
individual (GJ876) 

1. Introduction 

The search for habitable planets is essential in the quests to unravel the origin of the 
Solar System and find life elsewhere. Extra solar planets with mass {Mp) down to that 
of Neptune (Marcy et al. 2005; Santos et al. 2004) are thought to be gas giants because 
some have densities comparable to that of Jupiter. Extrapolation on the fraction of stars 
which host the yet-to-be-detected terrestrial planets depends on the assumed models 
of planet formation. If giants were formed through gas accretion onto pre-assembled cores 
(Pollack et al 1996), they would be accompanied by many earth-like planetary siblings. How- 
ever, in systems where giants were formed through GI (Boss 1997), such an association would 
be coincidental. 

Both models were introduced and modified to account for the observed properties of 
giants and extra solar planets. Several solutions (Hubickyj ct al. 2004; Alibcrt et al. 2005) 
have been suggested to bypass the protoplanctary growth bottlenecks associated with the 
early SA models. Likewise, photocvaporation has been adopted into the Gl model (Boss 2003) 
to address the giants' internal structure and correlation with metal-rich stars. While these 
revisions have led to a more comprehensive understanding of the planet-formation process, 
they have also reduced the capability to falsify the models. 

We present a robust observational test to differentiate the SA and GI models. In §2, we 
briefly recapitulate their essential features. Based on the SA model, we identify in §3 several 
potential avenues which may lead to the prolific formation of a population of close-in earths 
around stars with 1) short-period, 2) solitary eccentric gas giants, or 3) intermediate-period 
resonant planets. As a specific example, we suggest in §4 that the dynamical structure of the 
system (two giants with 30 and 60 day orbits plus a 7.5 M® planet with a 2-day orbit) around 
GJ 876 (with a mass = 0.32Mq) was established when the migrating resonant giants 
induced dynamical evolution of residual embryos and led to the formation of the close-in 
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earth. Since the embryos of this planet formed outside the snow hne, it is mainly made of 
water and it probably has an oxygen-rich atmosphere. In §5, we summarize the dominant 
SA processes which promote the formation of close-in earths. Since these planets are not 
expected to form if the giants originated through GI (Boss 1997), they provide an observable 
test for the SA and GI models. 



2. Differences between SA and GI models 

In the SA paradigm, planetcsimals coagulate via runaway growth and become proto- 
planetary embryos (with mass M^) by oligarchic growth (Kokubo & Ida 1998) on a time scale 
Tgrowth- III gaseous disks, the embryos' growth are stalled (Lissauer 1987) when they have 
swept clean the planetcsimals within their fccding-zone of full width Afz ~ 10 — 12Rh, where 
Rh = (Me/3M,)i/3a and a are their Hill's radius and semi major axis. For a fiducial surface 
density S^^g = So^_ot,/d,g(a/lAU)"^/^ [the subscripts d, g refer to the heavy elements and gas 
with scaling factors fa and fg with respect to the minimum mass nebular model (hereafter 
MMN) for the mass distribution inside 1 and 10 AU around the Sun (Hayashi et al. 1985)], 
the isolation mass is 




The normalizations, T^od = 107iccg cm ^ and = 2.4 x lO'^/depg cm ^ are determined 
by the volatile-ice enhancement [interior ['j^ce = 1) and exterior (7icc = 4) to the snow line 
(Ida & Lin 2004; Ida & Lin 2005) at ajce ^ 2.7{MjMQf AU] and the gas-depletion factors 
at time r* since the onset of star formation [fdep = exp(— r*/rdep) where the gas depletion 
time scale Tdep is observed (Hartmann 1998) to be a few Myr]. 

Prior to severe gas depletion, the embryos' eccentricities (ee's) are damped (Ward 1993) 
on a time scale re,damp - 300/g"Vd^p(^e/A^'e)~H«/lAU)2 yr << Tgrowth- Isolated embryos 
emerge on a time scale tiso = 3rgrowth, where 

„,o -1,-1,-2,". / « {'^'''A"'' f ^'•V"\, tn-, 

In disks with Tgrowth < Tdep, the embryos' growth at a < Oice is limited by Mjgo << M^. 
Just outside Oice, fa is enhanced by the grains' interaction with the disk gas (Stevenson & Lunine 1988), 
with Me > Mcore(~ 3 — lOM®) required for efficient gas accretion (Hubickyj et al. 2004) 
so that primary giants emerge there. However, in low-E,^ disks, no giants can form be- 
cause Tgrowth (-^core) > "^dep Gveu with a local fd enhancement. In modest-E,^ disks, when 
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Tgrowth(-^core) ~ ^dep, the rate of gas accretion and the asymptotic Mp are hmited by the 
global disk depletion or local gap formation (Bryden et al. 1999). Beyond the gap, grains 
and planetesimals accumulate which promote the emergence of secondary giants with a delay 
of Ar (Bryden et al. 2000). 

During the active disk-evolution phase, giants also migrate (Ida & Lin 2004) on a time 

scale 

Tmig - O.8/-Vde;(Mp/M^)(M0/M,)(lO-V«)(a/lAU)V2Myr, (3) 

where a is an ad hoc scaling parameter for angular momentum transfer efficiency. In persis- 
tent disks, extensive migration leads to short-period planets (Lin et al. 1996), but in rapidly 
depleting disks, they may be stalled. Multiple planets formed in slowly evolving disks 
(with Ar > Tmig) attain wide separations {eg Ups And). In rapidly evolving disks (with 
Ar < rjcp and < r^aig) the migration of successive emerging planets lead to resonant capture 
(Lee & Peale 2002) {eg GJ 876 and 55 Can). As they continue to migrate, the giants' ec- 
centricities (ep's) are excited by their resonances and damped by their tidal interaction with 
the gas beyond the gap (Kley et al. 2004) on a time scale Tp^damp ~ Tmig- 

In the SA hypothesis, giants form outside the snow line while many embryos with Migo 
remain in the inner disk. In the GI model, giants must form in massive disks which can effi- 
ciently cool (Gammie 2001; Durisen et al. 2005). These conditions are only satisfied in the 
optically-thin outer regions of disks around stars with r* < 0.1 Myr. Since the dynamical time 
scale of such regions and the growth time scale of the instabihty are ~ lO^yr, which is << 
Tgrowtii! the GI scenario requires terrestrial planets to form well after the giants (Boss 2003). 
In disks with fg » 1, the efficiency of angular momentum transfer is strongly enhanced 
by the growth of unstable modes (Laughlin & Rozyczka 1996; Armitage & Hansen 1999) so 
that a > 10~^ and rmig ~ 10^ yr. Giants' tidal interaction with gas in their co-orbital region 
of massive disks may also lead to rapid migration (Masset & Papaloizou 2003), though the 
magnitude of this effect remain uncertain (D'Angelo et al. 2005). Migrating giants would 
not capture each other into the observed mean-motion resonances (hereafter MMR) if their 
Tmig is shorter than their resonant libration time scale (Ida et al. 2000). 



3. Emergence of close-in ecirths 

With a series of simulations, we show that the embryos formed prior and interior to the 
giants are induced to migrate, collide, and evolve into close-in earths. We compute the evo- 
lution after the formation of the giants with a fourth-order time-symmetric Hermite scheme 
(Kokubo et al. 1998; Aarseth 2003) which includes the total force due to the star, planets, 
and disk (Mardling & Lin 2002; Nagasawa et al. 2003). The T,g distribution is assumed to 
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be that of the MMN with a gap (of width no less than the giants' Rh and epicycle ampli- 
tude). We adopt values of r^ig, Tdcpi, M^, and Mp to be within the range of the observed 
values (Hartmann 1998; Marcy et al. 2005). We highlight some generic features with three 
representative models: 

• a) a Jupiter-mass (Mj) giant formed at a = 5.2AU with Cp = (due to planet-disk 
tidal interaction (Papaloizou et al. 2001)) around a IMq star, 

• b) two giants [with initial (Mp, a, Cp) = (1.67Mj, 1.5AU, 0.40) and (3.1Mj, 4.17AU, 
0.24)] around a IMq star, 

• c) a pair of giants [with initial (0.56Mj, 3.3AU, 0) and (1.89Mj, 5.5AU, 0)] around a 
O.32M0 star, 

to represent systems around HD209458, /j, Ara, and GJ 876 respectively. 

In Model a), we adopt fg = 1 and fd = 1.25 which yield 27 embryos with initial 
Miso = 0.1 - 0.26Me and separation = lORh for a = 0.49 - 2.10 AU. (In all models, 
Miso ~ Mf-ore and Tgrowth(M:ore) < '^dep at o > Oice-) The orbits of cmbryos exterior to this 
range are rapidly destabilized by the giant whereas Miso of the close-in embryos does not 
contribute significantly to the total mass (Mtot = 4.35M0) of the population. 

The giant is assumed to migrate over Tmig — 1.8Myr (a = 10^"^) and to stall at 0.1 
AU. Embryos along the giant's migrating path are captured by its MMR's, migrate with it, 
and merge into four bodies with Mp — 0.11, 0.3, 1.09, & 1.92M0 just inside its asymptotic 
2:1 and 3:2 MMR's. Some captures also occur through the resonant interaction among the 
embryos. Since Cp = 0, this solitary giant's secular perturbation on the embryos is negligible. 
Due to the embryos' internal tidal dissipation (Novak et al. 2003) and relativistic precession 
(Mardling & Lin 2004), wc anticipate a slight additional ae decay and predict that short- 
period giants {eg HD 209458b) are accompanied by close-in earths (Figure 1). 

In Model b), we choose fg = fd = 4, and rdcp(~ O.lMyr < Tmig) with a = 10~^. The frac- 
tion of T Tauri stars with such massive disks (Beckwith & Sargent 1996; Muzerolle et al. 2005) 
is comparable to that of nearby stars with known planets. The initial 26 embryos with 
Miso = 0.23 - 1.25M® (a total of 14.7M®) are separated by = 12Rh at a = 0.09 - 1.2 
AU. The secular interaction between the two planets (with Mj and aj) introduces a preces- 
sion in their longitudes of periapse zup (initially orthogonal) with two eigenfrequencies which 
are also modified by the disk's potential (Ward 1993; Nagasawa et al. 2003). Due to similar 
effects and the star's post-Newtonian gravity, the embryos (with and mean motion rie) 
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also precess with the frequency 

^e-n^Y. ^aPf{a,)+2nU,Zik)^ + 3{^n (4) 

j=l,2 * ^ * 

Here aj — ae/aj, c is the speed of hght, and the power index of in MMN is A; = 3/2, for 
which the constant Z(k) — —0.54. Secular resonance (hereafter SR) occurs at locations 

where The resonant embryos' Cg are excited by the nearly aligned periapse (tUg — 

ZUp) and damped by the gas drag which also leads to orbital decay (Nagasawa et al. 2005). 
All precession frequencies and decrease with the disk depletion over Tdcp, leading to a 
sweeping SR. As their orbits cross, the embryos coagulate into two bodies with (Me,a)~ 
(7Me,0.23AU) and (4Me, 0.052AU) (Figure 2). 

In Model c), we adopt fg = fd = 1, Tdcp = 1 Myr, and Tmig = 1 Myr (a = 10^^) for the 
outer planet. The initial 24 embryos (with 0.04Me < < 1.49Me and a total ^ 7.2Me 
separated by = 12i?/i between 0.07-0.72 AU) extend outside the snow line (oice = 0.28 
AU). With an external disk, the outer planet undergoes orbital decay, captures the inner 
planet onto its MMR, and they continue to migrate together. The resonant excitation of the 
outer planet's Cp is damped by the outer disk (Kley et al. 2004) over Tpdamp ~ O.lr^i^. The 
inner planet's damping rate is assumed to be reduced by the gas depletion in the gap. 

Although the embryos' initial « ap, their Cg's are excited by the giants' pertur- 
bation and damped by their interaction with the disk. These effects lead to Oe decay 
(Nagasawa et al. 2005) on a time scale ~ e" ^re,damp > ^Vnig. As Op approaches their Oe's, 
embryos are captured into the gas giants' MMR's. The resonant embryos are induced to 1) 
migrate, 2) capture other embryos, and 3) coagulate over Tmig. After the giants are stalled 
at 0.2 and 0.12 AU, they continue to excite and reduce through their secular per- 
turbation and resonance, similar to Model b). After 1.4 Myr, three embryos remain with 
(Me,a)~ (1.5Me,0.072AU), (5Me, 0.054AU), and (0.7Me, 0.036AU) (Figure 3). Since we 
anticipate the giants' SR to sweep inward, embryos' eg-damping to weaken, migration to slow, 
and coalescence to occur, we continue the calculation with a reduced Tdep (0.1 Myr) and find 
two remaining embryos after 1.7 Myr with (6.5Me , 0.056AU) and (0.7Me , 0.037AU). The 
giants' secular perturbation on the embryos is suppressed by the relativistic precession which 
limits Cg ~ 0.01. The embryos undergo fractional orbital decay as they tidally interact with 
the host star during its life span (Mardling & Lin 2004). 

Extensive model analysis will be presented elsewhere. In general, the embryos' migration 
driven by the giants' MMR and SR arc robust (independently of fg, Tdep, 0(). However, the 
asymptotic values of Mg and Og arc determined by fd, Te,damp, Mp, the asymptotic values of 
Cp and ttp (which regulate Cg and embryos' response). 
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4. The origin of the close-in earth ciround GJ 876 

We choose Model c) to represent the initial disk structure around GJ 876. Both fg 
and fd in Model c) are high for disks around low- mass stars (Muzerolle et al. 2005). If 
Mdisk oc M^, {i.e. fd oc M^), Mjgo << and Tgrowth >> Tdep such that giants would 
rarely form around low-mass stars (Ida & Lin 2005). The total mass of refractory material 
interior to Oice is ~ fdM^. If the newly discovered planet is mostly composed of refractory 
material (~ 1.5% of all that in GJ 876), fd would need to be > 8 with a corresponding 
Mdisk > 0.02(/p//d)(a/lAU)^/^M* which is gravitational unstable at a > 1 - 2AU. 

Model c) reproduces the observed properties of GJ 876. The total amount of solids 
(mostly volatile ice) within lAU is ~ 7 — 8M^ and the condition for giant formation 
(Tgrowth(M;ore) < 7"dep 866 §2) is Satisfied at ~ 3AU (where the disk is gravitationally stable). 
This model suggests that the close-in earth was formed outside the snow line, mostly of 
water. While the detectable atmospheric H2O is photo-dissociated by the stellar UV flux 
(Leger et al. 2004) and hydrogen atoms escape, the atmospheric oxygen atoms are separated 
from the silicate core by a deep ocean. This process may lead to non-biogenic formation of 
O2, confusing the uniqueness of its bio-marker (Des Marais et al. 2002). 

5. Summciry 

The ubiquity of close-in earths is inferred from the SA hypothesis because their for- 
mation is promoted by: 1) the anticipated formation of embryos prior to the emergence of 
giants, 2) the driven migration of the embryos by the giants' MMR and SR, 3) the tidal inter- 
action between the embryos and their nascent disks, and 4) the embryos' induced collisions 
along their migration path. These effects are especially important for stars with multiple 
intermediate-period resonant planets. In the GI scenario, embryos cannot form prior to the 
giants. The giants' orbits evolve too rapidly for them to capture other giants and embryos 
onto their MMR's. Their detection through frequent radial velocity (Narayan et al. 2005) 
and sub-milli magnitude transit observations of known short-period and resonant giants 
could be used to extrapolate the probability of finding terrestrial planets around solar-type 
stars (Ida & Lin 2004). 

A general Hd.g distributions (other than MMN) may lead to diverse Miso, Tgrowth, and 
planetary configurations including earths with exterior to the Op of giants (Raymond et al. 2005). 
Provided T,d.g does not increase rapidly with a, close-in earths are expected to be associated 
with short-period giants. 

In Model c), aice's around M stars are much closer to those around G stars (in models 
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a and b). We infer from the SA scenario a prolific production of water- rich close-in earths 
around these low-mass stars. Although their day-side temperature may exceed 500K, the 
night-side of these close-in earths may be much cooler (Burkert et al. 2005) because their 
atmosphere is likely to be covered by opaque clouds. Direct searches for these conspicuous 
close-in earths may reveal that terrestrial planets are ubiquitous and lead to the detection 
of habitable environment among them. However, the possibility of a physical origin for O2 
requires a reassessment of the biological implications for its future detection. 

We thank G. Laughlin and S. Vogt for useful conversations. This work is supported 
by NASA (NAGS5-11779, NNG04G-191G), JPL (1228184), NSF (AST-9987417), NSFC 
(10233020 & NCET-04-0468), MEXT (MEXT 16077202), KITP & IPAM. 
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Fig. 1. — A solitary migrating giant captures residual embryos interior to the snow line onto 
its MMR. The time evolution of the semi major axes of the giant and embryos for Model 
a) is plotted with heavy and light lines. The dots and triangles represent collisions and 
final configurations, and the dotted fine shows the location of 2:1 resonance with the planet. 
Under the combined infiuence of ee-damping and the planetary resonant perturbation, the 
captured embryos migrate, collide, and grow. 
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Fig. 2. — Embryos' migration and collisions induced by the giants' SR in Model b). As SR's 
(dotted lines) sweep inward (Ward 1981) from Up and outward from the origin over r^ep , 
the migrating embryos collide and grow in agreement with theory. They stall when the SR 
is significantly weakened or is dominated by the relativistic correction. 
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Fig. 3. — The embryos' resonant capture and forced evolution by migrating resonant planets 
in Model c). As their orbits cross, the embryos undergo cohesive collisions and grow. The 
dotted line shows the location of 2:1 resonance with the inner planet. 



